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Abstract 



Possibility of establishment of equilibrium between the photon and the 

dense photon bunch is studied. In the case, when the density of plasma does 

not change, the condition of production of the Bose-Einstein condensate is 

obtained. It is shown that the inhomogeneity of density of photons leads to 

a new intermediate state of the photon gas. 
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Importance of scattering processes in interaction of photons with an electron gas was 
noted for the first time by Kompaneets J]]. Who has shown that the establishment of equi- 
librium between the photons and the electrons is possible through the Compton effect. The 
role of the Compton effect in setting equilibrium between the photons and the electrons was 
discussed by Kompaneets in the nonrelativistic approximation. In his consideration, since 
the free electron does not absorb and emit, but only scatters the photon, the total number of 
photons is conserved. Taking into account the fact that for the high frequencies of photons 
the probability of Compton scattering exceeds the absorption probability, Kompaneets has 
derived the kinetic equation for the photon distribution function. Using the kinetic equa- 
tion of Kompaneets, Zel'dovich and Levich |2| have shown that in the absence of absorption 
the photons undergo Bose-Einstein condensation. Such a possibility for the Bose-Einstein 
condensation occurs in the case, when the processes of change of energy and momentum in 
scattering dominates over the processes involving change of photon number in their emis- 
sion and absorption. Thus, in the above papers |[], [0 the authors demonstrated that the 
usual Compton scattering leads to the establishment of equilibrium state and Bose-Einstein 
condensation in the photon-plasma medium. 

In the present letter, we show that exists an another new mechanism of creation of 
equilibrium state and Bose-Einstein condensation in a nonideal dense photon gas. In our 
consideration we assume that the intensity of radiation (strong and super-strong laser pulse, 
non-equilibrium cosmic field radiation, etc.) is sufficiently large, so that the photon-photon 
interaction can become more likely than the photon-electron interaction. Under these con- 
ditions in the case of relativistic intensities of radiation of electromagnetic (EM) waves, we 
may consider the system of two weakly interacting subsystems: the photon gas and the 
plasma, which slowly exchange energy between each other. In other words, relaxation in the 
photon-plasma system is then a two-stage process: firstly the local equilibrium is established 
in each subsystem independently, corresponding to some average energies E 1 = KbT^ and 
E p = KsT p , where Kb is the Boltzmann constant, T 7 is the characteristic "temperature" of 
the photon gas, and T p is the plasma temperature. 
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As it was mentioned above, in the case of the strong radiation of EM waves there is an 
appreciable probability that the photon-photon interaction takes place, the phases of the 
waves will, in general, be random functions of time. We need therefore not be interested 
in the phases and can average over them. Under these conditions the perturbation state of 
the photon gas can be described in terms of the occupation number N(k, u, f, t) of photons 
and study how these numbers change due to the processes of interaction of the photons with 
others or with plasma electrons. Here, N(k,u,f,t) is the slowly varying function in space 
and time. 

A new version of kinetic equation for the occupation number N(k,u,r,t) of photons, 
for modes propagating with wavevector k and frequency u, at a position r and time t, was 
derived in Refs. IB, [1, B 



^N(k,u,r,t) + ^(k- V)N(k,u,r,t) - u 2 p sin^(v?- V fc " ;| " ^)/^ = > (!) 



where uj p = \/ 4 ^"° , f^oe is the electron rest mass, n and n are the non-equilibrium and 



equilibrium densities of the electrons, respectively, p = ^ ~, and 7 is the relativistic gamma 



factor of the electrons and can be expressed as 
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7~7>- 1 dk fduj N(k,u,r,t) 



ITiuj 

where 3 = ^ and % is the Planck constant divided by 2tt. 

Equation (|l]) is the generalization of the Wigner-Moyal equation for the classical EM 
field. Note that from this equation follows conservation of the total number of photons, but 
not the momentum and energy of photons, i.e., 

N = 2jdrJ —— J ^N(k, u, r, t) = canst. , (3) 

where coefficient 2 denotes two possible polarization of the photons. Hence, the chemical 
potential of the photon gas is not zero. 

In the geometric optics approximation (sinx ~ x) Eq.(^) reduces to the one-particle 
Liouville-Vlasov equation 



-iV^, W> r, *) + -(*■ V)N(k, u, f, t) - -f (v r -p • V ? - ^ • ^) 1 'J } = . (4) 



We specifically note here that in Eqs. (Q) and (§) there are two forces of distinctive nature, 
which can change the occupation number of photons. One force appears due to the redistri- 
bution of electrons in space, Vn e , and time, The other force arises by the variation of 
shape of wavepacket. In other words, this force originates from the alteration of the average 
kinetic energy of the electron oscillating in the rapidly varying field of the EM waves (V7, 



essentially changes on the distance L and in time t, which satisfy the inequalities u p t ^> 1 
and L ^> Under these conditions, the density of plasma almost does not change in 
comparison with the energy of EM field, i.e., we take into account the forces only due to the 
variation of 7. In this approximation, the statistical equilibrium between photons and the 
photon bunch will establish itself as a result of the " Compton" scattering processes. 

For convenience we introduce momentum and the "rest" energy of the photon as p = hk 
and m 7 c 2 = respectively. With these notations in the quasi-stationary case Eq.([|) 
takes the form 



The photon gas is characterized by an arbitrary initial distribution with some average 
energy. This may, in particular, be a Bose-Einstein distribution with a characteristic temper- 
ature T 7 , which, in general, will differ from the plasma temperature. Slower processes of the 
equalization of the photon and the plasma particle temperatures will take place afterwards. 

The solution of Eq.(||) can be written as 



dt 



In the following we shall consider the case, when the slowly varying function N(k, ou, r, t) 




(5) 





where /x 7 is the chemical potential of the photon gas. 

From Eq. (§) follows that the energy of the photon may be written as 
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e 7 = c^Jp 2 + m 2 c 2 . (7) 

The form of expression (^) coincides with the expression for the relativistic energy of material 
particles, so that the rest mass m 7 is associated with the photon in a plasma. We note that 
in contrast to material particles, the rest mass of photons depends on the plasma density and 
the relativistic gamma factor through the electron mass. Hence, the rest mass of photons is 
the function of the spatial coordinate r. 

The photon density n 7 , chemical potential, temperature and relativistic gamma factor 
are related by the expression 

where z = exp{ — m ~' c T ~ >1 ~' j is the fugacity, and = c^jp 2 + m 2 c 2 — m 7 c 2 is the kinetic 
energy of the photon. 

Requirement that the photon density must be positive leads to the important relation 
z e T i > 1, which should be satisfied for any energies, including the ground state, = 0, 
and at any points of space, i.e., the following inequality must hold 

m ^° 2 = > ^ ■ ( 9 ) 

The critical temperature of the Bose-Einstein condensation is determined for the fixed points, 
Tf, from the condition 

m 7(^/)c 2 = ^j^y = H • (10) 

After integration over the momentum, with condition /i 7 = m 7 c 2 , from Eq.([§D we obtain 
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n i = ^2^-Y. — K 2 (r}a) , (11) 

where A c = -^— c = ^-J 1 ^ 2 is the Compton wavelength of photon, a = K%{X) is the 

McDonald function of second order, and rj is an integer. 

We now first consider the non-relativistic "temperatures", i.e., when a ^> 1. Since, in 
this case K 2 {X) ~ \j§c) e ~ X > we ^ ne following expression for the critical temperature 
from Eq. fly]) 



2 1 / 3 7lh 2 2/3 

c ~ C 2 / 3 (3/2)m 7 ^ ' (12) 

where £(3/2) = 2.612 is the Rieman zeta function. 

For the ultrarelativistic photon gas, when the characteristic temperature satisfies the 
inequality m 0e c 2 ^> T 7 ^> m 7 c 2 , for the critical temperature we obtain 



T c = ^(^) 1/3 nf- (13) 



In the case, when the temperature T 7 is below the critical temperature T c of the Bose- 
Einstein condensation, the occupation number can be written in the form 

N(p, r) = + 4n 3 n ,6(p) . (14) 

ZXP{— t, J - 1 

Here, the second term corresponds to the finite probability for the photon to have exactly 
zero momentum, and no 7 is the density of photons in the ground state. 

Using Eq.([14]), we now determine density of the photons that fall into Bose-Einstein 
condensate. First, in the non-relativistic limit (m 7 c 2 > T 7 ), we have 

n ~f = n 7 (l - (jjff 1 ) , T 7 < T c . (15) 

Next, for the ultrarelativistic photon gas (m 0e c 2 > T 7 > m 7 c 2 ), we get 

n 07 = n 7 (l-(^) 3 ) , T 7 <T C . (16) 

It should be noted that the production of the Bose-Einstein condensates is realized, as follows 
from Eqs. (|9|) and fllOD, in regions, where the field intensity is large. 

Now we will show that the inhomogeneity of intensity of the EM field leads to the 
intermediate state of the Bose-Einstein condensates of the photon gas. The point is that the 
relativistic gamma factor, 7, is the function of spatial coordinates and at different points 7 
is different, which means that the effective mass m 7 of the photon, as a function of 7, also 
changes in the usual space. Hence, one can find some points or regions of high intensities 
in the gas volume, where equality (0) is fulfilled, so that the photon gas is found in the 



Bose-Einstein condensate state, whereas in the rest of the domains or points the photons 
remain in normal state, for which inequality (|9|) is satisfied. We can therefore conclude that 
in such a case the photon gas, as whole is in the intermediate state. This effect resembles 
the Landau's intermediate state in superconductor. 

In order to make this statement more clear, we consider periodically modulated cylindri- 
cally shaped laser pulse (light bullet). We assume that the field intensity is uniform in the 
transverse direction (in general, it is not necessary) of propagation of pulses and periodic 
along the light bullet (ID optical lattices). In this case the relativistic gamma factor has 
the following form 



where Qq is a constant, d and ^ are the depth and the period of modulation, respectively. 

Furthermore, we assume the periodic boundary conditions N(p,z) = N(p,z + L), and 
7(2) = 7(2 + L), with L being the length of the light bullet. 

We then suppose that at all points z — —I (where I = 0,1, 2, 3, ...) photons are in the 



Bose-Einstein condensate state, that is /x 7 = TiiOpj Jl + Qq(1 + d). Obviously, at other points 



Thus, the light bullet is divided into numerous parallel alternating thin layers of the 
normal and the Bose-Einstein condensate photon gas. This state of the photon gas may be 
called the intermediate state. 

We should emphasize that this phenomenon can also exist in the case of the usual Comp- 
ton processes, when the photon-electron interaction is more important than the photon- 
photon one. In this case, in the expression (|5|) we should take into account the variation of 
density of the electrons. 

In conclusion, we have presented a new concept of establishment of the equilibrium 
between the photon and the wavepacket of the EM field (the dense photon bunch). This is 
the fully relativistic effect due to the strong EM radiation. We have studied the possibility 
of creation of the Bose-Einstein condensates for the case, when the density of plasma does 




(17) 




the condition ( [TDD i s n °t valid, so that the photons stay ordinary. 
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not change. We have shown that the inhomogeneous dense photon gas can be found in 
the intermediate state. We think that this state is most common in interaction of the 
relativistically intense EM waves with plasmas. We believe that the intermediate state can 
be experimentally observed in atomic gases of alkaline elements contained in magnetic traps 
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